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ABSTRACT 

We have observed 28 local galaxies in the wavelength range between 1 and 2.4 



> 

' fj,m in order to define template spectra of the normal galaxies along the Hubble se- 

I quence. Five galaxies per morphological type were observed in most cases, and the 

. resulting RMS spread of the normalized spectra of each class, including both intrinsic 

' differences and observational uncertainties, is about 1% in K, 2% in H and 3% in J. 

Many absorption features can be accurately measured. The target galaxies and the 
spectroscopic aperture (7"x53") were chosen to be similar to those used by Kinney et 
al. (1996) to define template UV and optical spectra. The two data sets are matched 
Q^i in order to build representative spectra between 0.1 and 2.4 ^m. The continuum shape 

I ' of the optical spectra and the relative normalization of the near-IR ones were set to 

O I fit the average effective colours of the galaxies of the various Hubble classes. The re- 

sulting spectra are used to compute the k-corrections of the normal galaxies in the 
^ ' near-IR bands and to check the predictions of various spectral synthesis models: while 

the shape of the continuum is generally well predicted, large discrepancies are found 
in the absorption lines. Among the other possible applications, here we also show how 
these spectra can be used to place constraints on the dominant stellar population in 
local galaxies. Spectra and k-corrections are publicly available and can be downloaded 



' from the web site http://www. arcetri. astro. it/'^filippo/spectra. 
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1 INTRODUCTION 

Galaxies are complex systems where gas, dust, stars and 
dark matter deeply interact with each other for a very long 
time. The emerging spectrum is therefore sensitive to a 
number of parameters as mass, age, metallicity, star for- 
mation history, dust, geometry etc. Despite this complexity, 
many groups have tried to model the galaxy emission, and 
the observed optical spectra can now be accurately repro- 
duced (e.g. Bruzual & Chariot, 1993; Worthey, 1994; Fioc & 
Rocca-Volmerange, 1997; Barbaro & Poggianti, 1997; Tan- 
talo et al., 1996). This is increasingly important as more 
high-redshift galaxies are discovered and important quanti- 
ties as age, mass and star formation history are estimated 



by comparing the integrated properties of the galaxies with 
the predictions of spectrophotometric models. 

These models are calibrated on local galaxies and the 
first test on their properties is the comparison of the mod- 
eled spectra with the observed ones over a large wavelength 
range. Until now, this comparison could only be done in 
the optical and near-UV region where representative galaxy 
spectra are available (e.g.. Pence, 1976; Coleman, Wu & 
Weedman, 1980; Kennicutt, 1992a, 1992b; Kinney et al., 
1996, hereafter K96), while photometric points are used for 
all the other wavelengths. The aim of this work is to extend 
the range where a full comparison is possible towards longer 
wavelengths, between 1 and 2.4 fim. This is in particular 
important because large differences are present between the 
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different models especially in this wavelength range where 
the relative contribution of dwarf, giant and supergiant stars 
is poorly known (e.g., Chariot, Worthey & Bressan 1996). 
These stellar populations are characterized both by differ- 
ent depths of the absorption hues and different shapes of the 
near-IR continuum (see, for example, Pickles, 1998; Langon 
et al., 1999) and therefore by studying breaks, colours and 
lines of a full spectrum it is possible, in principle, to separate 
their contribution. 

Near-IR spectra can also be used to compute the k- 
corrections of an unevolving galaxy population, and, there- 
fore, to study the luminosity functions and the differen- 
tial evolution of the galaxies (Pence, 1976; Coleman et al. 
1980, Frei & Gunn 1994, Poggianti, 1997). The near-IR 
k-corrections are becoming more and more important as 
new IR-optimized telescopes are becoming available, as the 
Very Large Telescope, the Large Binocular Telescope and 
the Next Generation Space Telescope. Complete spectra be- 
tween UV and near-IR wavelength can also play a role in the 
measure of the redshifts of large samples of galaxies, both 
from photometric data (for a review see Yee 1998) and even 
from spectroscopic data (eg., Glazebrook, Offer and Deeley, 
1998). 

In this paper we show the results of the observations of 
28 galaxies of morphological type between E and Sc, while 
later classes will be the subject of a future work. Target 
selection, observations and data reduction are described in 
sec. ^. Great care was taken in matching the near-IR spec- 
tra to the optical ones by using the observed UV-optical-IR 
colours, as described in sec. ^ and sec. ^. In sec. |^ we derive 
the k-corrections in the J, H and K band. In the following 
sections, we compare the total spectra to the predictions of 
some spectrophotometric models for elliptical galaxies and 
use the observed absorption lines to study the dominant 
stellar populations. Spectra and k-corrections are available 
in electronic form. 



2 OBSERVATIONS AND DATA REDUCTION 

The target galaxies were selected in the morphological 
classes E, SO, Sa, Sb and Sc. In each class we selected large 
nearby galaxies without indication for peculiar activity. An 
absolute magnitude threshold of Mv < —21 was used to 
ensure a homogeneous metallicity and to allow a sample as 
close as possible to that in K96. We also mainly observed 
galaxies whose optical spectra were available in the literature 
to be also able to study single objects. Therefore we selected 
galaxies from the samples by K96 and Kennicutt (1992a, 
1992b), adding one more elliptical galaxy from the list given 
in Goudfrooij (1994). The resulting 28 targets are listed 
in Table |l], together with the morphological informations 
derived either from the LEDA catalog ( jittp : / /leda. univ^ 
lyonl.fr/leda/| ) or from the reference listed in the table. The 
objects are sorted in the table and divided into 5 classes 
according to the morphological index T (see, for example, 
Buta et al., 1994). 

We planned to have 5 galaxies for each class between 
E and Sc with good S/N ratio (SNR) to have an estimate 
of the intrinsic spread in the spectra of each class and to 
average out peculiar characteristics. This number was not 
reached for the Sa class because NGC2798 turned out to 



Table 1. Observed galaxies 



NGC 


R.A. (1950) DEC. 


Type 




Ref.'' 


Bands'^ 


E (-5 <T<-3) 


4648 


12:39:54.5 


4-74:41:44 


E3 


-4.9 


1 


JHK 


1700 


04:54:28.1 


-04:56:30 


E4 


-4.8 


2 


JHK 


3379 


10:45:11.3 


-1-12:50:48 


E2 


-4.8 


1 


JHK 


4472 


12:27:13.9 


-1-08:16:22 


El 


-4.7 


1 


JHK 


4889 


12:57:43.7 


-1-28:14:54 


E4 


-4.3 


1 


JHK 


SO (-3 <T<0) 


1023 


02:37:15.5 


-1-38:50:56 


E-SO 


-2.6 


3 


JHK 


3245 


10:24:30.1 


-1-28:45:45 


SO 


-2.1 


1 


JHK 


4350 


12:21:25.1 


-1-16:58:21 


SO 


-1.8 


3 


JHK 


4382 


12:22:52.8 


-1-18:27:59 


SO-a 


-1.3 


3 


JHK 


5866 


15:05:07.0 


-1-55:57:20 


SO-a 


-1.3 


1 


JHK 


Sa (0<T<2) 


2681 


08:49:58.0 


-f51:30:14 


SO-a 


0.4 


3 


JHK 


3623 


11:16:18.6 


-1-13:22:00 


SBa 


1.0 


1 


JHK 


2775 


09:07:41.0 


-1-07:14:35 


Sab 


1.7 


1 


JHK 


3368 


10:44:06.9 


-1-12:05:05 


SBab 


1.7 


1 


JHK 


Sb (2<T<4) 


4826 


12:54:16.9 


-1-21:57:18 


Sab 


2.4 


3 


JHK 


4736 


12:48:31.9 


-1-41:23:32 


Sb 


2.5 


3 


JHK 


2841 


09:18:34.9 


-1-51:11:19 


Sb 


3.0 


3 


HK 


4102 


12:03:51.6 


4-52:59:23 


SBb 


3.1 


3 


JHK 


3147 


10:12:39.3 


4-73:39:02 


Sbc 


3.7 


1 


JH 


Sc (4<T<6) 


2903 


09:29:20.2 


4-21:43:19 


SBbc 


4.0 


1 


JHK 


5194 


13:27:45.9 


4-47:27:12 


Sbc 


4.1 


3 


JHK 


3994 


11:55:02.3 


4-32:33:23 


Sc 


4.9 


3 


JHK 


1637 


04:38:57.5 


-02:57:11 


Sc 


5.0 


3 


HK 


2276 


07:10:22.0 


4-85:50:58 


Sc 


5.4 


1 


K 


Others 


3432 


10:49:42.9 


4-36:53:09 


Irr 


9.0 


3 




2798 


09:14:09.4 


4-42:12:34 


Sap 


1.7 


1 




1569 


04:26:05.8 


-1-64:44:18 


Irr 


8.2 


1 




4449 


12:25:45.9 


4-44:22:16 


Irr 


9.3 


1 





": morphology stage index (e.g., Buta et al., 1994) 

**: references. 1: Kennicutt (1992b); 2: Goudfroij (1994); 3: K96. 

'^: bands used for the templates 



be peculiar (it is part of an interacting system) and with 
emission lines. This galaxy is listed in the "Others" section 
of Table ^ together with 3 more observed objects of later 
stages not used for the templates. 

All the spectra were obtained at TIRGO, an IR- 
optimized 1.5m telescope on the Swiss Alps. We used the 
long-slit spectrometer LonGSp (Vanzi et al., 1997b) which 
is based on a NICMOS3 256x256 pixels Rockwell array and 
is sensitive between 0.9 and 2.5 /xm. 

Each galaxy was observed through the largest available 
slit, 7" wide, in order to have a good coverage of the galaxy. 
Spectra were obtained in the J, H and K bands for inte- 
gration times between 15 and 30 minutes for each band. 
The resulting spectral resolution is about 300 in J, 400 in 
H and 500 in K. The telescope was chopped every minute 
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1.5 



G 



0.5 



1.15 1.2 1.25 1.3 

Figure 1. Rest frame average spectra of each class of galaxies in 
the J band. The thick line is the average of the observed spec- 
tra, the thin lines show the ranges within 1 standard deviation. 
Arbitrary offsets were added to the spectra for clarity. 



between the galaxy and the nearby sky to have a good sky 
subtraction. Each time the galaxy was placed on a different 
position of the slit to minimize the effect of bad pixels and 
distortions of the flat-field. For each galaxy, stars of spectral 
type between F8 and G8 were observed in order to mear 
sure the atmospheric transmission and the variation of the 
instrumental efficiency with the wavelength, as discussed in 
Maiolino & Rieke (1996). The total number of acquired im- 
ages is about 7000. 

All the spectra were sky subtracted using the average 
of the adjacent spectra, and divided by differential flat-fields 
obtained by imaging the dome. The resulting frames were 
then rectified, co-aligned and stacked together using a clip- 
ping algorithm. A one-dimensional spectrum was then ex- 
tracted from the central 53" in order to maximize the ob- 
served fraction of the galaxies. The resulting aperture of 
7"x53" is similar to the aperture of 10"x20" used by K96 
to define their templates and allows us to merge the two data 
sets. At the average distance of the galaxies of our sample, 
Mpc (^50 =Ho/50 km/sec/Mpc), the spectroscopic 
aperture corresponds to about 1x7 Kpc. We estimated the 
fraction of the galaxy light collected by the slit by using the 
existing aperture photometry in the V band (Prugniel and 
Heraudeau, 1998): on average, our spectra contain about 
20% of the galaxy light, ranging from 7% for the closest and 
largest galaxies to about 35% for the smallest ones. There- 
fore our results should not be considered "total" spectra, 
even if the difference is expected to be negligible in the near- 
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A (/im) 

Figure 2. As Figure 1, H band 
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Figure 3. As Figure 1, K band 
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Figure 4. K-band spectra of the galaxies with active star forma- 
tion NGC1569, NGC2798 and NGC4449, The emission lines at 
2.122/im (H2) and 2.166 (Br7) can be seen in the spectra. 



IR (see discussion below) . This is more important when also 
the K96 spectra are considered as in the optical the radial 
gradients are expected to be larger (but see also the discus- 
sion in K96). 

The "true" spectra of the reference stars were modeled 
by the solar spectrum (Livingston & Wallace, 1991) cor- 
rected for the slightly different temperature of the stars. The 
derived transmission curve was then applied to the galaxy 
spectra. The final spectrum of each galaxy in each band is 
then corrected for Galactic extinction, normalized, shifted 
to zero redshift, and averaged with the other galaxies of 
the same class. A clipping algorithm was applied to remove 
a few residual discrepant value due to cosmic rays, erratic 
bad pixels not present in the mask and uncorrected features 
of the reference stars. In Figures |l|,^ and ^ we show the aver- 
age spectra and the relative standard deviations in the three 
bands. 

The differences between the spectra of the galaxies of 
the same class is due to both all the observational uncertain- 
ties (noise, residual sky contribution, distortion on the flat 
field, non perfect correction of the instrumental and atmo- 
spheric transmission) and the intrinsic differences between 
the observed galaxies. We can therefore have a robust upper 
limit of the total uncertainties of our spectra by looking at 
the resulting RMS among the galaxies of the same type. This 
quantity is between 0.7% and 1.4% in K, between 1.7% and 
2.4% in H and between 2.4% and 3.6% in J, depending on 
the morphological class. The tendency to have larger spreads 
towards bluer wavelengths is probably dominated by larger 
intrinsic differences, and the value in K, about 1%, can be 
assumed as the real uncertainty in the single spectrum. Be- 
cause of the different redshift of the galaxies, less galaxies 
contribute to the edges of the spectra and therefore the un- 



certainties tend to be larger in the external part (about 5%) 
of the bands. The parts of the spectra without indication of 
standard deviation in figures 1,2 and 3 have been derived 
from only a galaxy and should be treated with caution. 

If the normal galaxies show a high degree of uniformity, 
on the contrary large differences are found for the few galax- 
ies with peculiar properties or of later types. The K spectra 
of 3 of these galaxies (NGC 1569, NGC2798 and NGC4449) 
are shown in Figure while for NGC3432 the SNR is too 
low. The galaxies in the figure show prominent emission lines 
whose brightness and ratios change much from one object 
to the other. In this case no attempt was made to create a 
template. 



3 CALIBRATION OF THE SPECTRA 

We now want to establish the relative calibration of the spec- 
tra in the various bands to be able to study the overall spec- 
tral energy distribution from the UV to the near-IR. It is 
not possible to derive these relative normalizations by sim- 
ply using spectral observations because of several effects, 
such as slit losses, variability of the atmospheric transmis- 
sion and uncertainties on the reference star colours. The 
usual method to normalize spectra taken in different bands 
is to use the broad-band colours of the target objects: spec- 
tra should be calibrated to reproduce the colours of the ob- 
served galaxies inside the spectroscopic aperture. 

As this work is aimed at the construction of spectra rep- 
resentative of the average properties of the different classes 
of galaxies, we have calibrated the spectra by using the av- 
erage colours of each class. This method guarantees that the 
final spectra closely match the average properties of the class 
and allows the reduction of the uncertainties on the colours 
as hundreds of galaxies can be used. 

Of course, colours inside our rectangular spectroscopic 
aperture are not available. However, the dependence of near- 
IR colours and spectra with aperture is known to be very 
weak because the dominant stellar population at these wave- 
length is always the giant branch (with the exception of star- 
burst galaxies, see for example, Maraston (1998)) and only 
the second-order metallicity and age gradients could be seen. 
Therefore we can use colours inside circular apertures as 
they are the only ones usually published, and the difference 
in shape between photometric and spectroscopic apertures is 
negligible if the two areas are chosen to be similar. The most 
similar circular aperture with enough published photometric 
data is the "effective" one, i.e., the aperture containing half 
of the light of the galaxy in the B band. Effective colours are 
available for a large number of galaxies and in many filters 
and is therefore possible to have large homogeneous data 
sets. 

For the smallest objects, the effective aperture is similar 
to the spectroscopic one, for the largest galaxies the photo- 
metric aperture is larger. We have performed two checks to 
be sure that, even in the latter case, effective colours can reli- 
ably be used to calibrate the spectra: first we have compared 
spectra extracted from the central 10" of the slit with those 
of the outer part; second, we have compared the spectra of 
large and small galaxies of the same class (as NGC4350 hav- 
ing Do=2'.0 and NGC4382 with Do=8'.2) . In both cases no 
systematic differences are seen. 
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Table 2. Average effective colours of galaxies with <-21 
For each colour the standard deviation and the number of used 
objects are also reported. 



U-B 


E 


0.50 




(0.08) 




323 


SO 


0.47 




(0.11) 




287 


Sa 


0.36 




(0.19) 




138 


Sb 


0.22 




(0.20) 




321 


Sc 


0.06 




(0.18) 




294 


Sd'' 


-0.12 




(0.16) 




53 




-0.15 



B-V V-R 



V-I 



V-K J-H" H-K" 



(0.20) 
102 



0.99 
(0.05) 
418 

0.97 

(0.08) 
344 

0.90 
(0.11) 
185 

0.82 
(0.12) 

541 

0.70 
(0.13) 
536 

0.62 
(0.18) 
99 

0.51 
(0.17) 
117 



0.59 
(0.05) 
314 

0.58 
(0.05) 
227 

0.58 
(0.08) 
73 

0.57 
(0.09) 
156 

0.52 
(0.10) 
133 

0.47 
(0.13) 
25 

0.40 
(0.20) 
28 



1.22 
(0.07) 
221 

1.20 
(0.08) 
158 

1.17 
(0.11) 
82 

1.16 
(0.11) 
315 

1.15 
(0.15) 
287 

1.09 
(0.19) 

58 

1.08 
(0.30) 

35 



3.30 
(0.09) 
32 

3.25 
(0.14) 
13 

3.24 
(0.18) 
17 

3.21 
(0.28) 
16 

3.03 
(0.24) 
23 

2.95 
(0.32) 
12 

2.35 
(0.35) 
5 



0.66 
(0.05) 
225 



0.21 

(0.02) 
225 



0.66 0.22 

(0.05) (0.02) 

235 235 

0.67 0.25 

(0.06) (0.03) 

105 105 



0.66 
(0.06) 
93 

0.66 
(0.07) 
46 



0.25 
(0.03) 

93 

0.25 
(0.04) 
46 



0.65 0.23 
(0.08) (0.05) 
26 24 



0.51 
(0.10) 
22 



0.21 
(0.06) 
20 



°: The J-H and H-K colours are based also on the results in Fioc & 
Rocca-Volmerange, 1999, where only average quantities are given. 
In these cases the scatter is not measured but estimated. 
Mv <-20 

no magnitude selection. 



Once the aperture is chosen, we have also to select the 
luminosity of the galaxies used to derive average colours. 
As described in the previous section, the target objects 
are large, metal rich galaxies, and we want to use galaxies 
with the same properties to calibrate their spectra. In fact, 
smaller galaxies tend to be less metallic and therefore bluer 
in the optical colours (see, e.g., Fioc & Rocca-Volmerange, 
1999), giving origin to the well-known colour-magnitude re- 
lation. We will use optical and near-IR colours appropriate 
for galaxies with Mv <-21, corresponding to our sample. 

Many photometry papers where used to derive these 
colours and check the results (Aaronson, 1978; Frogel et 
al., 1978; Persson et al., 1979; Griersmith et al., 1982; 
Glass, 1984; Giovanardi & Hunt, 1988; de Jong & van der 
Kruit, 1984; de Vaucouleurs et al., 1991; Buta et al., 1994, 
1995a, 1995b; Prugniel & Heraudeau, 1998, Fioc & Rocca- 
Volmerange, 1999). The optical colours are mainly based on 
the Prugniel & Heraudeau (1998) catalog after correcting for 
Galactic dust using the RC3 value of the galactic extinction 
in B (de Vaucouleurs et al., 1991) and the Cardelli et al. 
(1989) extinction curve. The E and SO optical-to-IR colours 
are based on Persson et al. (1979) and Frogel el al. (1978), 
while the colours of later spirals are mainly derived from 
Griersmith et al. (1982), de Jong & van der Kruit (1984) 




- 



Figure 5. Colour correction to the K96 spectra. The dotted lines 
show the K96 spectrum for the ellipticals and a spline fit to the 
corresponding photometric points, (empty dots). The solid dots 
show the effective colours of the galaxies of this class with My < 
—21 (see Table 1^), the solid lines are the fit to these data and the 
final spectrum, corrected to fit the colours. 



and Aaronson (1978). IR-to-IR colours are based on Glass 
(1984), Frogel et al. (1978) and Fioc & Rocca-Volmerange 
(1999). Note that the J-H and H-K colours of the early-spiral 
galaxies are redder than those of the ellipticals: a short dis- 
cussion about this unexpected effect can be found in Fioc & 
Rocca-Volmerange (1999). All the colours were converted to 
Johnson's U, B and V, Cousins' R and I and CIT J, H and 
K. Great care was used in comparing different galaxy sam- 
ples: whenever possible, the results from different samples 
were compared to each other to check for consistency and to 
discover any important selection effect. The results are listed 
in Table ^ together with the colour scatter and the number 
of galaxies used. The absolute magnitude threshold was re- 
duced to M\/=-20 for the Sd class (6<T<8) because these 
galaxies are on average intrinsically fainter than the other 
classes, and for the I class (T>8) no magnitude restriction 
was used because of the large scatter observed. The colours 
of these latter two classes (not used for the templates and 
reported for completeness) are based on a smaller number 
of galaxies and are probably less accurate that the others. 

Earlier classes tend to be more homogeneous than the 
later ones. For the ellipticals, the spread of the colour dis- 
tribution, when removing a few very discrepant objects, is 
comparable to the uncertainties on the single measure as 
quoted in the original papers, leaving little room to an in- 
trinsic spread of the population due, for example, to diflerent 
metallicity or star formation histories. On the contrary, for 
the later class the observed spread is several times larger 
than the measured uncertainties. As an example, for the 
ellipticals the V-K colour distribution has a standard devi- 
ation of 0.09 mag (and an error on the mean of the order 
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of 1%), similar to the average estimated error on the single 
point of about 0.10 mag, while for Sc the observed spread is 
0.24 mag. 



4 MATCHING OPTICAL AND NEAR-IR 
SPECTRA 

The same colour calibration was used to match our spectra 
to existing optical and UV ones to obtain a single template 
from 0.1 to 2.4 /xm. Several optical templates are available 
in the literature (Coleman et al., 1980; Kennicutt, 1992a, 
1992b; K96). In order to obtain meaningful results, the opti- 
cal spectra must have appropriate resolution and wavelength 
coverage, should be based on a relatively large number of 
galaxies and the used spectroscopic aperture must be simi- 
lar to our one. Among the available optical templates, only 
the K96 spectra have the required features: 1- they have a 
high enough resolution (between 6 and lOA) for most studies 
of stellar populations; 2- have the largest wavelength cover- 
age, extending from 1200 A up to about 1/^m, at the edge 
of the J band (the Sc spectrum covers only up to 7700 A); 
3- use an aperture of 10x20" at all the wavelengths, similar 
to ours. 4- are based on the observations of 30 galaxies, 12 
of which in common to our sample. 

While below 7000A the K96 spectra closely reproduce 
the observed colours of Table ^ above this wavelength the 
K96 spectra are much redder. The difference is not due to the 
chosen colours (effective aperture and < —21) because 
the K96 spectra don't fit the average colours inside any aper- 
ture and for any magnitude (see, for example, Goudfrooij, 
1994). Moreover, they don't reproduce the colours of the ob- 
served galaxies: as an example, all 4 elliptical galaxies in K96 
have observed Johnson's V-I colours between 1.38 and 1.45 
for apertures between 10 and 20" (Prugniel & Heraudeau, 
1998), while the colour derived from the K96 spectrum is 
1.59. It should also be noted that above 7000A the K96 spec- 
tra are quite different from those by Coleman et al. (1980) 
who also corrected the spectra for the observed colours. No 
useful comparison can be done with the Kennicutt (1992a, 
1992b) atlas because it only covers the spectra region below 
7000A. 

In order to have a self-consistent template we have to 
apply some corrections to the optical spectra to make them 
fit the colours. The procedure is shown in Fig. ^ where the 
K96 spectrum for the ellipticals is plotted against the colours 
of Table ^ used both to calibrate the near-IR spectra and to 
correct the K96 templates. The correction to the K96 spectra 
was computed by the ratio of two splines, fitted respectively 
to the observed colours and to the colours derived from the 
spectra (see Figure ^). This correction is not unique, but the 
availability of data in many filters maintains this uncertainty 
below the other ones. 

For all the K96 spectra between E and Sb, a good agree- 
ment (within a few percent) is present below 7000A between 
the optical spectra and the average colours, and the correc- 
tions are small. On the contrary, the Sc spectrum, based on 
two galaxies only, appears to be much bluer than the ex- 
pected colours of its class even at A < 7000 A. Also the 
equivalent widths of the emission lines are much larger than 
in the galaxies of our sample, about 296A for the Ha-|-[NII] 
vs. an average of 4lA for our galaxies (Kennicutt & Kent, 



1983) and 54A for the Sc sample by Kennicutt (1992a). We 
think that both these effects are probably due to the fact 
that the two Sc galaxies in the K96 sample, NGC 598 and 
NGC 2403, are of later type (T=5.6 and 5.7) then the av- 
erage Sc (T=5.0, from Buta et al., 1994) and those in our 
sample (T=4.7). To correct for this, we have both applied 
the colour correction, as for the other types, and also scaled 
down the emission line to have the same EW of our sam- 
ple. The resulting spectrum, although representative of an 
average Sc galaxy, should be used with caution when com- 
paring optical to near-IR features and when dealing with 
the emission lines, also because of its lower SNR. 

As final step, the spectra were arbitrarily interpolated 
across the wavelength ranges where the atmosphere is not 
transparent. In Figure ^ we show at a reduced resolution the 
final spectra reproducing the observed effective colours. 

The relatively small apertures used for both spec- 
troscopy and photometry means that the spectra in Figure ^ 
are representative of the central regions of the galaxies con- 
taining about half of the total emission. K96 compare their 
optical spectra with those by Kennicutt (1992b) obtained 
using a much larger aperture (about 90x90") and conclude 
that for E, SO and Sa the K96 spectra are in fact repre- 
sentative of the total emission of the galaxies because their 
integrated properties are dominated by the central bulge. 
The same argument applies in the near-IR too, and in this 
case a weaker dependence of the spectra from the aperture 
is also predicted: in fact, while in the optical the young pop- 
ulation in the upper main sequence can easily dominate, 
in the near-IR any shortly living super-giants population 
must compete with the usually much more numerous giants 
population and therefore dominates only in case of strong 
activity of star formation. As discussed in the sec. ^ this 
can be checked comparing near-IR spectra extracted from 
the inner and outer regions of the galaxies, or from large 
and small objects. A weaker variation of the near-IR spec- 
tra with aperture with respect to the optical spectra can also 
be predicted from the colours: the J-H and H-K colours are 
only weakly dependent of the galaxy type (see Table |^ and 
do not show significant radial gradients (Glass 1984; Fioc & 
Rocca-Volmerange, 1999), not even for the late type galax- 
ies. Therefore we conclude that for the early type galaxies 
the spectra in Figure |^ can also be used as templates of the 
integrated galaxy. For Sb and Sc galaxy, younger extra nu- 
clear populations can change the spectra significantly, both 
in the optical because of the stars of the upper main se- 
quence and in the near-IR if an important population of red 
supergiants is present. 



5 THE K-CORRECTIONS 

The first quantities that can be derived from these spectra 
are k- correct ions, i.e., how magnitudes and colours change 
with the redshift because of the shift of the observed wave- 
length ranges. These corrections are crucial to link the prop- 
erties of the local universe with the distant one (see, for ex- 
ample, Peebles, 1993), and their accurate knowledge is im- 
portant to study the evolutionary effects. Our spectra allow 
us to compute these corrections in the near-IR filters start- 
ing from z=0 with a greater accuracy than by using simple 
broad-band photometry as, for example, Cowie et al. (1994). 
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Figure 6. Template spectra for E, SO, Sa, Sb and Sc galaxies. The spectra were reduced to a resolution of 40A and normalized to the 
B band. The interpolation between the bands is arbitrary. 



In Figure ^ we show the k-corrections computed for the J, 
H and K bands. These filters have an increasing importance 
in the surveys for distant galaxies and in the study of the 
known high-redshift objects, especially when more 8-meters 
class space- and ground-based telescopes optimized for the 
infrared will be available. 

The k-corrections were computed using Bessell & Brett 
(1988) filters, more standard that the CIT system. We note 
that in the K band the k-corrections are almost independent 
from the galaxy type up to z=2, reflecting the dominance of 
giant stars in these wavelengths in all the observed galaxies, 
as previously noted (Glazebrook et al. 1995). Note that no 
contributions are included from galaxy evolution or increas- 
ing absorption from the intergalactic medium; the computa- 
tion is limited to z=2 because above this limit these effects 
are certainly dominant. 

We have compared our k-correction in the K band with 
others previously published. At low redshift {z < 0.6) we 
found consistency with the k-correction by Glazebrook et 
al. (1995) derived from Bruzual & Chariot (1993) model for 
a SSP with age of 5 Gyrs, and used, for example, by Loveday 
(2000) and Szokoly et al. (1998) to define the local luminos- 



ity function. On the contrary, big differences (up to 0.2 mag) 
are found with k-correction from 'UV-hot' elliptical model of 
Rocca-Volmerange & Guiderdoni (1988). At higher redshift 
{z > 0.6) our k-corrections are 'redder' up to 0.3 magnitudes 
than previous ones, however at these redshifts they should 
be taken with caution, and considered as lower limits, due 
to the possible evolutionary effect on galaxy spectrum. 

Important differences are also found with the k- 
corrections by Poggianti (1997). Regarding the SED of the 
elliptical, the difference in J is generally about 0.1 and in 
H is always below 0.2, while differences greater than 0.3 are 
found in the K band in the redshift range z=0.3 to 1.3, up to 
a value of about 0.45 magnitudes at z=0.4-0.7. Although the 
differences in the filter response functions adopted here and 
in Poggianti (1997) contribute to the differences between the 
K-corrections, most of the K-band discrepancy arises from 
intrinsic differences in the SED, probably due to the sparse 
sampling of the Poggianti (1997) models at A > 17000 A. 
The comparison for Sa and Sc galaxies is less meaningful 
because of large spectral variations within each morpholog- 
ical class, possible aperture effects (see §4), especially for Sc 
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Figure 7. k-corrections in the J, H and K bands between rcdshifts 
and 2 for each morphological type: E (solid line), SO (dotted), 
Sa (short dashed), Sb (long dashed) and Sc (dot-dashed). Our 
spectra define these k-eorrections up to z=0.3 in J, z=0.7 in H 
and 1.2 in K, while at higher redshift they reflect the K96 spectra 
corrected to match the observed colours. 



galaxies, and the presence of emission lines, not included in 
the Poggianti (1997) model. 



6 COMPARING OBSERVATIONS AND 
MODELS 

Template spectra can be used to test the spectrophotomet- 
ric models of the galaxies. Despite their low resolution, the 
SNR is high enough to allow useful comparison of both the 
continuum shape and the spectral features. Here we show 
two examples of the possible use of these spectra. The first 
one, described in this section, is the comparison of the shape 
of the continuum over a large wavelength range. In the next 
section the absorption lines will be compared. 

Two reasons suggest to start the continuum compari- 
son from the elliptical galaxies: first, from the observational 
point of view this class is the most homogeneous, its aver- 
age colours are better defined and the spectrum has smaller 
uncertainties. Second, elliptical galaxies arc known to be 
dominated by old stellar populations and therefore the un- 
certainties on the star formation history (SFH) have lower 
consequences than in the case of the spirals: while in the 
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Figure 8. Mean deviation of the predictions of some spectropho- 
tometric models from the observed spectrum of the elliptical 
galaxies as a function of the age of the stellar population. Details 
about the models can be found in the text. 



latter case detailed SFH arc required to fit the spectra, for 
the ellipticals a single parameter, the age, is usually enough 
to obtain a good agreement. 

We used four different models: the Bruzual & Char- 
lot model in the 1999 version (BC99), the Fioc & Rocca- 
Volmerangc (1997) PEGASE model (FRV97) the Barbaro 
& Poggianti (1997) one (BP97), and the model by Worthey 
(1994) (W94). We have also used a previous (1996) version 
of the Bruzual & Chariot model (BC96) that allows for non- 
solar metallicities. 

BC99 and BC96 use isochroncs from the Padova group 
(Bressan et al., 1993) and two different stellar libraries, the 
empirical library from Pickles (1999) and the theoretical one 
by Lejeune et al. (1997). The latter library permits to pre- 
dict the spectra of galaxies of various metallicities, from 5% 
to 250% solar. BC99 and BC96 give the spectra of an instan- 
taneous burst of star formation (the so called Simple Stellar 
Population, SSP) that can be integrated over the time to re- 
produce an arbitrary SFH. We used a Salpeter Initial Mass 
Function (IMF) between 0.1 and 125 M©. The BC96 model 
was used mainly because it allows to use non solar metallic- 
ities. 

The FRV97 models use a different approach: starting 

from some hypothesis on the efficiency of the star formation 
processes, they model the entire life of the galaxy in order to 
reproduce the observed spectrum. The resulting SFH for el- 
lipticals is a rapidly declining exponential law giving most of 
the star formation in the first Gyr, while for the Sc the SFH 
is almost constant over an Hubble time. The chemical evolu- 
tion is also computed, but is not used for the final spectrum 
which assumes solar metallicity. The Padova isochrones are 
used, together with the empirical stellar library by Gunn & 
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Figure 9. The observed spectrum of the elliptical galaxies (thick line) is compared with a SSP spectrum of the BC99 model (thin line) 

with an age of 12 Gyr and solar mctallicity. The upper panel shows the full spectrum, the four lower panels show enlargements in the 
optical and in the three ncar-lR bands. The shape of the spectrum is very well reproduced in all the range between 0.2 and 2.4 fj,m, and 
this is true for I lie oirtical absorption features too. In the near-IR large discrepancies can be noted, especially in J and H. The reason of 
these differences are explained in the text. 



Stryker (1983) and Langon & Rocca-Volmerange (1992) and 
the IMF by Rana (1992). 

The BP97 model has the same approach of BC96, com- 
puting the galaxy integrated spectrum adding up the con- 
tribution of SSPs of (possibly) different metallicities. It uses 
the Padova isochrones and the Salpeter IMF between 0.1 
and 100 M^. The stellar spectra in the optical arc cither 
from Kurucz (1993) or from Jocoby et al. (1994), while the 
library by Langon & RoccarVolmerange (1996) are used in 
the near-IR. In this paper we only consider BP97 spectra 



up to 1.8 /um because above this limit the Kurucz sampling 
becomes too sparse. 

The W94 model provides with single-age stellar popu- 
lation with metallicities between 0.01 and 3 times solar. It 
uses stellar evolution isochrones by VandenBerg and collab- 
orators and by the Revised Yale Isochrones (VandenBerg 
& Laskaridcs, 1987; Green et al., 1987), Salpeter IMF and 
theoretical stellar spectra by Kurucz (1993) and BcsscU et 
al. (1994). A particular emphasis was put in studying the 
variation with metallicity and age of several optical spectral 
features. 
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For demonstration purposes and to better study the 
intrinsic differences between the models, for BC99, BC96, 
W94 and BP97 we used SSP spectra of solar metallicity. 
For FRV we used both the models of elliptical galaxy they 
present, named E13 and E16, corresponding to two different 
cosmologies. In all the cases we searched for the best fitting 
spectrum by changing the age, the only remaining parame- 
ter, between 5 and 19 Gyr. All the spectra, both observed 
and modeled, were normalized to have the same total flux 
between 0.3 and 2.35 /^m and reduced to the same resolu- 
tion. The difference between comparing spectra to spectra 
and spectra to broad-band colours is important: both the 
absorption features and the continuum shape inside each 
band are taken into account and can dominate the results, 
especially when the agreement with the colours is good. 

To measure the deviation it is not possible to use the 
standard because the fltting points are not totally inde- 
pendent of each other as the normalization of each band is 
given by the colour calibration (see section ^) . Therefore we 
need a way to compare the deviations both to the spectrum 
noise and to the colour uncertainties. To do this we compute 
the square deviation A'^ between observed {oi) and model 
(mi) spectra in each band: 



1=1 



■mi) 



N 



where A'^ is the number of spectral points, and divide it for 
the square sum 5 of the various error contributions. The 
final value A we use as analog to is the average of these 
ratios in the optical and in three IR bands: 
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which has a value of 1 if the deviations are only due to 
the errors in the observed spectra. The quantities 5^ are 
difflcult to estimate with precision: from the RMS spread of 
the spectra (see section ^ and the error on the mean of the 
colours in table ^, we obtain 5 = 0.03 in the optical, J and 
H and 5 = 0.02 in K. 

In Figure ^ we show the resulting values of yl as a func- 
tion of the galaxy age for the four used models and solar 
metallicity. Even taking into account the uncertainties on 5, 
for all the models an age exists which well reproduce the 
data: BC99 shows the best accord for ages between 10 and 
16 Gyr, W94 select younger ages, between 8 and 9 Gyr, 
FRV97 reaches the best agreement for older ages, above 16 
Gyr. The deviation of the BP97 model is always above 1 but 
steadily reduces toward older ages, reaching a value of about 
2 at 16 Gyr, the maximum age available for the model. It 
should be noted that BC99, BP97 and W94 are single-burst 
models with the same IMF and metallicity, therefore the 
difference between them are due to something outside the 
user control, as isochrones or stellar spectra. FRV97 is not a 
single burst model, about 92% of the stars are formed in the 
first Gyr (about 97% in the first 4 Gyr) but the SFR never 
reaches a null value; this probably explains at least part of 
the differences with the other models and the requirement 
of very old ages. 

In Figure ^ we show the best fltting model, the BC99 
spectrum for an SSP of 12 Gyr of age, solar metallicity and 
empirical star spectra. The overall shape of the spectrum is 



very well reproduced: the largest difference are in the inter- 
polated region between 1 and l.l/im where the model trace 
the position of the continuum while our interpolation an "ef- 
fective" level taking into account the absorption lines. On 
the contrary, at a closed look it is possible to see that impor- 
tant discrepancies remains in the details of the absorption 
lines, both if theoretical and observed star spectra are used. 
This will be the subject of the next section. 

We also tested the effect of changing the metallicity by 
using the BC96, W94 and BP97 models with metaUicities 2.5 
times solar, BC96 0.2 solar and W94 0.3 solar. In all cases no 
good fit can be reached for any age, as the values of A remain 
above 5. Models with super-solar metaUicities are too red for 
A > 8000 A, even for young ages (limited to be higher then 
5 Gyr): as an example, all these models predict fluxes in the 
H band between 10% and 20% higher than observed. Vice 
versa for sub-solar metaUicities: the H band flux is under 
predicted of about 20%. These discrepancies might be an 
indication of the real metallicity of the elliptical galaxies, 
but it is more probably due to inadequacy of the non-solar 
stellar library or isochrones. It is beyond the purpose of this 
paper to solve this ambiguity. 



7 THE ABSORPTION LINES 

The use of our spectra to study the absorption line is ham- 
pered by the low resolution, both because faint lines are 
below the detection threshold and because many lines are 
blended together. Nevertheless, many features can be eas- 
ily seen and measured with great accuracy, even if in some 
cases their identification may not be unique. For the fol- 
lowing analysis we have only used the spectra in H and K 
because a series of facts make our J spectra less useful to 
study the features: they have the lowest resolution, their 
SNR is lower than that in the other bands, the features are 
intrinsically less prominent, and the star classification in this 
band is still in its infancy (see, for example, Wallace et al., 
2000). 

To further increase the SNR we decided to use the av- 
erage of the spectra of E, SO and Sa because they show 
no significant difference in the absorption lines, meaning 
that temperature, luminosity and metallicity of the domi- 
nant stellar populations are similar. We will indicate with 
ES the resulting average spectrum. 

We identified the lines (see Table ^ by comparing the 
features in our spectra to those identified in the late type 
stars. The literature of reference is now rich enough (e.g., 
Origlia et al., 1993; Oliva et al., 1995, 1998; Meyer et al, 
1998, hereafter M98; Kleinmann & Hall, 1986, hereafter 
KH86; Vanzi & Rieke 1997; Wallace et al, 1996, 1997) to 
make the identifications reliable, but the dominant contri- 
bution to some features due to blending of many lines is 
sometime uncertain. In these cases, the synthetic spectra 
constructed specifically for line identification in the H band 
by Origlia et al. (1993) were of particular help, while in the 
K band the high-resolution spectra by Wallace et al. (1996) 
were taken into particular account. 

The first useful results can be derived by qualitatively 
comparing the H- and K-band parts of the ES spectrum to 
various stellar spectra libraries and to the results of the spec- 
trophotometric models. This comparison is shown in Fig- 
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Table 3. Feature Identification 



A 


Main 


Other 


(Mm) 
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1.529 


OH 


CN,TiI 


1.540 


OH 


Sil 


1.558 


i^CO 


OH 


1.577 


12CO 


MgI,FeI 
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Sil 


OH 


1.598 


12CO 


SiI,"CO 


1.606 


OH 




1.619 


i^CO 


OH,CaI 


1.640 


12CO 


SiI,[FeH] 


1.652 


"CO 


OH 


1.661 


12CO 


OH 


1.669 


Sil 


OH 


1.672 


All 


HI,i2cO 


1.677 


All 




1.689 


OH 


HI,CO 


1.710 


Mgl 


CO,OH 


1.723 


OH 


Sil 


1.733 


Sil 


HI 


2.067 


Fel 




2.072 


Pel 




2.081 


Fel 


Sil 


2.107 


Mgl 


H20,SiI 


2.117 


All 


H2,MgI,FeI 


2.136 


Sil 




2.146 


Mgl 


NaI,SiI,CaII 


2.166 


Br7 


VI 


2.173 


Scl 


Fel 


2.179 


Til 


SiI,FeI 


2.189 


Sil 


TiI,FeI 


2.208 


Nal 


ScI,TiI,VI,FeI,SiI 


2.226 


Pel 


ScI.Til 


2.239 


Fel 


Scl 


2.248 


Fel 


VI,TiI 


2.263 


Cal 


ScI,TiI,FeI,SI 


2.281 


Mgl 


CaI,FeI,SI,HF 


2.294 


i^CO 


Til 


2.323 


12CO 




2.345 


"CO 





ures and0 where all the spectra were normalized to a 
fitted continuum and reduced to the same resolution of the 
observed data. 

In panel (a) of Figure ^ a K5III star from M98 is over- 
plotted to the data. The agreement is very close and basi- 
cally all the features are well reproduced: as expected, cold 
giant stars dominate the emission at 1.6 /xm and the ob- 
served spectra can be reproduced even with a single star 
spectrum Analogous excellent agreement is obtained using 
spectra from the Wallance & Hinkle (1996) library (not 
shown) which are limited to a narrower wavelength range. 

The following panel shows the spectrum of two stars, 
a K5III and a MOIII, from the more recent stellar libraries 
by Pickles (1998), library used in several spectrophotomet- 
ric models. The former star, the same as panel (a), shows a 
spectrum similar to our data, even if the agreement is not 
as accurate as for the KH86 spectrum. On the contrary, the 
MOIII star shows a very different behavior, as also other stars 
of similar spectroscopic and luminosity class. These differ- 
ences are not explained by the small temperature difference 



between the K5 and MO stars (about 160 °K, M98) and such 
a behavior is not seen in other libraries, as M98 or Wallance 
& Hinkle (1996). The differences are probably due to the 
use of several different sources in the derivation of the Pick- 
les (1998) library. In particular Pickles bases the near-IR 
part of its spectral on the the Langon & Rocca-Volmerange 
(1996) library because it covers a wide wavelength range, 
between 1.43 and 2.5 /^m, with accurate flux calibration. 
Unfortunately, some features of the spectra of this library 
seem very variable from one spectral type to the next ones 
and very different from other libraries (as M98 and KH86) 
and from our spectra. Therefore they probably are observa- 
tional artifacts. As a consequence, any spectrophotometric 
model based on this library will not be able to accurately 
reproduce the features in near-IR galaxy spectra. This is 
shown in the following (c) panel: here the observed spec- 
trum is compared to the BC99 SSP galaxy model at the age 
of the best fit in Figure ^ This model is the only one in our 
sample having a resolution in the near-IR high enough to al- 
low a meaningful comparison of the spectral features. Even 
if the general shape of the continuum is well reproduced, 
as shown in Figure ^, the agreement with the spectra fea- 
tures is generally poor: as an example, the deep absorption 
lines predicted by the models at 1.665 and 1.705 ^m are not 
present in the observed data. 

Same conclusions follow from the spectra in the K band 
shown in Figure |l^: the K5III star from KH86 in panel (a) 
shows a very good agreement with the ES spectrum. Only 
the CO bands at 2.29 and 2.32 fim are deeper in the star 
spectrum than in the ES one, because of either metallicity 
effects, or the contribution of stars warmer than K5 or on the 
main sequence, or differences in the fitted continuum level 
above 2.28 ^m. Analogous excellent agreement is obtained 
with the K5III star spectrum from the Wallance & Hinkle 
(1996) library (not shown). 

In panel (b) and (c) the Pickles (1998) and BC99 spec- 
tra are shown: again, some features are not reproduced. For 
example, the prominent Nal feature at 2.21 fim is almost 
completely absent from the MOIII star and from the model 
spectra. Also the Cal feature at 2.26 /im appears deeper in 
the observed data than in the model spectrum. Given that 
the general shape of the continuum is well reproduced by 
the model, (see the previous section) the problem must be 
in the used stellar library. 

To summarize, libraries of stellar spectra exist which 
reproduce the observed galaxy spectra, while others show 
very different features. Spectrophotometric models based on 
the second set of spectra might not be able to reproduce the 
observed features. Our spectra will make it possible to tune 
the models, especially when new and more extended libraries 
of near-IR stellar spectra will be available (see, for example, 
Ivanov et al, 1999; Mouhcine & Langon, 1999). 



8 THE DOMINANT STELLAR POPULATION 

As discussed in the previous section, indications on the dom- 
inant stellar population can be easily derived even by simply 
comparing the observed data with the stellar spectra. More 
quantitative results can be obtained by comparing the ab- 
sorption lines. 

Several authors have used this method to classify the 
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Figure 10. Line identifications and comparison of the observed 
H-band spectrum of the early-type galaxies (thick line) normal- 
ized to its continuum to several spectral libraries (thin lines) re- 
duced to the same resolution. Panel (a): the galaxy spectrum is 
compared to a K5-giant star from the libraries by M98, showing 
a very good agreement. Panel (b): the spectra of a K5 (solid) and 
a MO giant (dotted) from Pickles (1998) are overplotted to the 
galaxy data. Even if the two stars are similar, their spectra in this 
library are very different. Panel (c): the SSP model from BC99 at 
10 Gyr with empirical stellar libraries (see previous section and 
Figure Bj) is compared to the data. 



stars by their near-IR spectra (see references below). Usu- 
ally two sets of lines are considered to derive both the spec- 
tral and the luminosity class, the first set depending both on 
temperature and gravity, while the second one on tempera- 
ture only. The recipes used for star classification can also be 
applied to the observed galaxy spectra to derive indications 
on the the dominant stellar population. 

We took into account two schemes of classification. 
They are both based on the CO(2,0) band head at 2.29 
Hm whose intensity depends both on temperature and lumi- 
nosity and can be used to study the stellar populations in 
galaxies or in dusty regions (e.g., Armus et al., 1995; James 
& Seigar, 1999; Mobasher & James, 2000). The two meth- 
ods use different thermometers: the first one, introduced by 
KH86 and refined by Ramirez et al. (1997), is based on the 
Equivalent Width (EW) of Nal 2.21/im and Cal 2.26^im 
which depend strongly on temperature but not on luminos- 
ity. In the second one, based on Origlia et al. (1993), the 
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Figure 11. As the previous figure, but in the K band. The K5III 
star spectrum in panel (a) is from KH86. 

temperature is measured by the ratio of the two features at 
1.62 and 1.59 /j,m. In both cases, the results for the galaxies 
are compared with those for the stars of known properties. 

To apply these methods, we measured the EW of the 
lines with respect to a local continuum obtained by fitting 
a straight line to the "clear" parts of the spectrum. The 
spectral region used for both the line and the continuum 
are listed in Table ^. The same procedure was used to mea- 
sure the EWs of the stars of the M98 and KH86. Because 
of the lower resolution of our spectra, our definitions are a 
little different from those by, for example, KH86 or Origlia 
et al. (1993), and the measured values of the EWs are con- 
sequently different. The errors on these quantities were esti- 
mated by taking into account both the noise of the spectra 
and the uncertainties in the continuum level. The latter con- 
tribution is usually the dominant one, especially in H where 
many features are present and it is difficult to define the 
continuum level. The results for the lines used in the clas- 
sification are shown in Table ^ for the average ES spectra 
and in Table ^ for the individual galaxies not used for the 
templates. It should be noted that the error on the 1.62/1.59 
ratio shown in the table is smaller than the error on each 
single line because the same continuum fit is used: the uncer- 
tainties on the continuum level strongly affect the measure 
of the EW of each line, but partially cancel out from their 
ratio. 

In Figure O we present the classification based on the 
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Table 4. Feature definition 



Feature 


Line 


Continuum 




(/^m ) 


{^J,m) 


1.59 


1.586-1.594 


1.572-1.574, 1.608-1.613 






1.628-1.633 


1.62 


1.616-1.629 


1.572-1.574, 1.608-1.613 






1.628-1.633 


Nal 


2.204-2.211 


2.191-2.197, 2.213-2.217 


Cal 


2.258-2.269 


2.245-2.256, 2.270-2.272 


12CO(2,0) 


2.289-2.302 


2.270-2.272, 2.275-2.278 






2.282-2.286, 2.288-2.289 



Table 5. Value of the EW (A) measured on the templates, and 
their errors 



Feature 


E-fSO-fSa 


Sb 


Sc 


1.59 


3.9 (0.8) 


3.7 (0.9) 


3.3 (0.9) 


1.62 


6.0 (0.7) 


5.6 (0.8) 


5.7 (0.8) 


1.59/1.62 


1.5 (0.1) 


1.5 (0.2) 


1.7 (0.2) 


12CO(2,0) 


14.6 (1.7) 


15.8 (1.7) 


15.6 (1.7) 


Cal 


2.6 (1.1) 


2.7 (1.1) 


2.2 (1.1) 


Nal 


3.4 (0.7) 


4.2 (0.7) 


3.4 (0.7) 



features in the K band studied by KH86 and Ramirez et 
al. (1997). This figure differs from both works because the 
feature definition is different (because of the lower resolu- 
tion) and because we don't take into account the EW of Br7 
which cannot be reliably measured. As expected, the results 
show that in all the cases the main role is played by class 
III stars, possibly with some contribution from dwarf stars. 
With the present data it does not seem possible to measure 
this contribution and therefore constraint the ratio between 
dwarf and giant stars, but this could be done with higher 
resolution observations. The spectra of the active galaxies 
(see Figure ^ are too noisy to measure the EWs of Nal and 
Cal, and only upper limits can be derived. As shown in Fig- 
ure ^ these upper limits put two of these galaxies in the 
region between class I and class III stars, suggesting that 
the contribution supergiants might be important in these 
galaxies. 

In Figure ^ we apply the second method based on the 
CO feature at 2.29/.im and on the ratio of the 1.59 and 
1.62/im in the H band. This method is based on the work 
by Origlia et al (1993) who showed that this ratio can be 
use as a thermometer in case of giants and supergiants stars. 
Also in this case the line ratio is consistent with the class 

Table 6. Value of the EW (A) measured for the active galaxies, 
and their errors 



NGC2798 


NGC1569 


NGC4449 


12CO(2,0) 11.8 (2.4) 


16.7 (2.3) 


13.8 (2.9) 


Nal+Cal < 3.7 


< 4.2 


< 5.8 




EW["^CO(3,0)] (A) 



Figure 12. The sum of the EWs of Nal (2.21 ^tm) and Cal (2.26 
fim) is plotted versus the EW of the 1^00(2,0) band at 2.29 fim. 
The solid dots with the error bars are our galaxies, while dwarf, 
giant and supergiant stars are plotted as empty squares, stars and 
empty circle respectively (see text for references). The lines are 
linear fit to the stellar data, with the indication of the luminosity 
class. 

Ill stars. 

A discussion about the metallicity of ellipticals and 
early spirals is beyond the scope of this paper. Detail discus- 
sions about the information of metallicity and element ratios 
based on optical data can be found, for example, in Worthey 
(1998) and Trager et al. (2000). Origlia et al. (1997) have 
shown that moderate-resolution near-IR spectra can also be 
used to measure the metallicity of the systems. Their recipe 
is based on the 1.62/im feature which is shown to depend on 
the metallicity index [Fe/H] and, more weakly, on the carbon 
depletion factor [C/Fe]. It also depends on the microturbo- 
lent velocity ^ of the dominant stellar populations because 
this quantity influences the strength of the saturated molec- 
ular lines. By observing globular clusters of known metallic- 
ity, Origlia et al. (1997) show that this metallicity scale gives 
errors smaller that 0.3 dex. These authors have also observed 
a limited number of elliptical galaxies. The few metallicity 
estimates published for these objects are based on the Mg2 
index (see, for example, Davies et al., 1987) and on the Mg2- 
metallicity calibration by Casuso et al. (1996). While the 
Mg2-derived value of [Fe/H] of the observed galaxies was 
solar or slightly supersolar, between -1-0.08 and -1-0.21, the 
values based on the near-IR lines turned out to be subso- 
lar, between -0.28 and -0.46. This can be explained either 
by the presence of a corresponding Magnesium enhancement 
of about [Mg/Fe]~0.5, as proposed by several authors (see 
Origlia et al., 1997 for details), or by an anticorrelation be- 
tween [C/Fe] and [Fe/H], the more metallic objects being 
more carbon depleted, as observed in the globular clusters. 
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Figure 13. EW of the ^^00(2,0) feature at 2.29/xm vs. the ratio 
EW(1.62)/EW(1.59). Giants and supergiants stars (see text for 
references) are plotted by empty circles and stars respectively. 
The galaxies observed by Origlia et al. (1993) are plotted by solid 
diamonds if quiescent ellipticals and spirals, and solid squares if 
HII galaxies. Solid circles with error bars lare the galaxies of this 
study. The lines are the linear fit to the stellar data. 

More recently, several authors (Trager et al., 2000; Worthey, 
1998) would rather explain the non-solar [Mg/Fe] with a de- 
pletion of the iron-peak elements. 

Despite the large uncertainties on the right value of 
[C/Fe] and on the resulting metallicity, we have applied 
this method to our early-type galaxies (ES spectrum, see 
table By iteratively solving the eq. 2 by Origlia et al. 
(1997) we obtain a mean microturbolent velocity ^ of about 
3.1 km/sec, in agreement with the values between 2.3 and 
3.0 in Origlia et al. (1997). By using the 'standard' value 
of [C/Fe] =— 0.3 suggested by these authors and derived for 
the Magellanic Clouds clusters by Oliva & Origlia (1998), 
the metallicity derived from their eq. lb is [Fe/H]~ —0.20, 
reducing to [Fe/H]~ -0.4 for [C/Fe]=0.0. The error on this 
quantity is difficult to estimate and is probably dominated 
by the uncertainties on [C/Fe], but this result can be used 
as an indication of sub-solar abundance of iron in the early- 
type galaxies. 



9 SUMMARY 

The near-IR spectra of 28 nearby galaxies were used to de- 
fine the "average" J, H and K spectra of normal galaxies 
along the Hubble diagram between E and Sc. The target 
objects and the aperture were chosen to be similar to those 
by K96 in order to merge the two data sets and obtain self- 
consistent spectra between 0.1 and 2.4 /jm. The average ef- 
fective colours of the galaxies of the various classes were used 



to calibrate the final spectra. The overall precision is about 
2%. 

The spectra were used to compute the k-corrections in 
the near-IR bands and to check the predictions of several 
spectrophotometric models of galaxy evolution. We find that 
most model can accurately reproduce the observed contin- 
uum shape of the elliptical galaxies, if solar metallicities are 
assumed. In this case the derived SSP ages at above 8 Gyr. 
Nevertheless, large differences are found between the mod- 
els, for example in terms of the derived age of the dominant 
stellar population. 

On the other hand, larger discrepancies are found be- 
tween the observed spectral features and the predictions of 
the models. We note that the spectra of cool giant stars from 
some accurate libraries (as KH86 and M98) can reproduce 
the observed features in great detail, and conclude that the 
disagreement is probably due to the stellar libraries used in 
the models. 

The spectral features are also used to study the stel- 
lar populations by using several stellar classification recipes 
based on near-IR absorption lines, some of them applied to 
galaxies for the first time. The dominant contribution to the 
near-IR emission of all the galaxies between E and Sc is con- 
firmed to be due to giant stars of spectral type about K5, 
with solar or slightly subsolar metallicity. 

The final spectra for each galaxy type, the relative k- 
corrections (both in the Bessel & Brett and in the CIT sys- 
tems) and the profiles of the used filters can be downloaded 
from the web site www.arcetri.astro.it/~filippo/spectra. In 
some cases it could be useful to simulate the spectra of galax- 
ies with colours different from those in Table |^, for exam- 
ple to reproduce individual galaxies or match total colours. 
In the same web site it is possible to define a new set of 
colours, compute the corresponding spectra and download 
the results. 
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